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Abstract 

We calculate the process-independent fragmentation functions for a b antiquark to 
fragment into longitudinally and transversely polarized B* ( 3 S\) mesons to leading 
order in the QCD strong coupling constant. In the special case of equal quark mass 
we recover previous results for the fragmentation of c — > ip and b — ► T. Various 
spin asymmetry parameters are defined as measures of the relative population of the 
longitudinally and transversely polarized vector meson states. In the heavy quark 
mass limit my, — > oo our polarized fragmentation functions obey heavy quark spin 
symmetry, we therefore apply them as a model to describe the fragmentation of charm 
and bottom into heavy-light mesons like D* and B* . The spin asymmetry parameter, 
a(z), is consistent with the existing CLEO data for D*. The scaling behavior of (z) 
is studied in detail. We find excellent agreement between the predictions of (z) from 
our fragmentation functions and the experimental data for D* and B* from the LEP, 
CLEO, and ARGUS detectors. Finally, we also point out that the spin asymmetry 
depends significantly on the transverse momentum p± of the vector mesons relative to 
the fragmentation axis. 
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I. Introduction 



The physics of hadrons containing a single heavy quark has been studied intensively in 
the past several years, mainly due to the development of the powerful technique of Heavy 
Quark Effective Theory (HQET) [|. O ne crucial observation in HQET is that the 

heavy quark spin decouples from the strong interaction dynamics in the limit of infinite 
heavy quark mass. This happens because the leading operator, the chromo-magnetic dipole 
moment, that couples the heavy quark spin to the gluon field is inversely proportional to the 
mass of the heavy quark. 

Whenever a heavy quark is produced with high degree of polarization (for example, 
bottom and charm quarks are produced with 94% and 67% left-handed polarization, respec- 
tively, at the Z° pole), it might be possible to extract the spin information of the heavy 
quark if the subsequent hadronization does not lead to substantial depolarization. The en- 
ergy spectrum or the angular distributions of the lepton in the semi-leptonic decays of the 
top, bottom, or charm quarks are often used as the spin analyzer of the heavy quark spin. 
Several recent works (Refs.Q, [5], |7], ||) have been devoted to the polarization effects in the 
D* and D** mesons, B* and B** mesons, and ££, and A& baryons at the LEP energy. In 
this article, we will investigate in detail the spin alignment of a heavy quark that undergoes 
fragmentation into the spin triplet S-wave meson state. 

Recently, it has been pointed out || that the dominant production of (be) mesons at the 
large transverse momentum region is due to fragmentation, in which a high energy b quark 
is produced from a hard process and subsequently fragments into various (be) meson states. 
Furthermore, the corresponding process-independent fragmentation functions of b — ► B^Sq) 
and b — > B*( 3 Si) were calculated within perturbative QCD (PQCD) || [Kj to leading order 
in both the strong coupling constant a s and v, where v is the typical velocity of the charm 
quark inside the meson. A lower bound of the inclusive branching fraction for the production 
of B c at Z° has been estimated to be about 2.3 x 10~ 4 0, including both IS and 2S states. 
It implies that only about 230 B+ or 230 B~ are produced from 10 6 Z°. At the Tevatron 
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with a luminosity of 25 pb -1 , one expects [11] about 2 x 10 4 B c mesons to be produced with 
Pt{B c ) > 10 GeV from the direct b quark and induced gluon fragmentation. The three- 
charged leptons from a secondary vertex observed in the decay — > J/ip + t'v$t followed 
by J/ip — > U, where £, £' = e, /i, can provide a clean signature for B c . The combined 
branching ratio of these decays is expected to be ~ 0.2%, which implies 40 distinct events 
at the Tevatron. In addition, the B c meson can be fully reconstructed via hadronic decay 
modes, e.g., B+ — > J/ijj + tt + , with J/xjj — > ££. Thus, unless LEP can increase its luminosity 
by an order of magnitude or so in the near future, the best place to look for the B c meson 
will be at the Tevatron. 



In Ref. [12 1, the asymmetry due to the relative probabilities for the production of the 
transversely versus longitudinally polarized J / ip states by c and c quark fragmentation was 
discussed. As a result of this asymmetry, an anisotropic angular distribution of the leptons 
in the decay J/ip — > ii was found to be of order 5%, which might not be large enough to 
be observed due to the presence of another important source of highly polarized J / ip from 
.B-meson decays. Since the heavy quark fragmentation probabilities, first calculated in [|TH 



are known to have big cancellation in the equal mass quarkonium case, one expects to have 
larger asymmetry in the unequal mass case like the B* meson system. Such asymmetry 
can be measured in principle from the anisotropy of the photon angular distribution in the 
decay B* —>■ B c + 7. To predict the degree of anisotropy one needs to know the polarized 
fragmentation functions for the b quark splitting into various helicity states of the B* meson. 

In this work, we will study the spin alignment in the B* production via the fragmentation 
of the b quark. First, we will calculate the fragmentation functions for b — > polarized B* 
states, as an extension to the previous calculation of b — > unpolarized B* M, or as an 
unequal mass extension to the calculation of c (or c) — > polarized J/ip |12| . This extension 
is also useful beyond the (be) system as the final results of the fragmentation functions can 
be applied phenomenologically to the case of heavy-light mesons such as D* and B*. This 
allows us to get a better insight to the spin asymmetry in the D* and B* mesons, where 
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perturbative methods are not applicable. 

Spin alignment of the B*, J/ip, Y, D*, and B* systems will be studied. Two spin asym- 
metry parameters frequently quoted in the literature are 

and a= — - — , (1) 



> L + T T 

where T(L) denotes the production probability of the transverse (longitudinal) state of the 
excited S-wave meson (B*, J/ip,Y,D*, or B*). Here we introduce another parameter 

T 

W = — . (2) 

T + 2L K } 

Two useful relations of these spin asymmetry parameters are the following 

, «=^. (3) 

Note that we did not specify the production mechanism in the definitions of these spin 
asymmetry parameters. 

To measure these spin asymmetry parameters, one can study the two body decay of the 
excited meson, e.g. B* — > Bj, D* — > Dtt, and B* — > B c ^ etc. The angular distribution of 
the emitted photon (or pion) depends on the helicity of the parent meson. For definiteness 
we will consider the B* meson in the following. Suppose the B* meson is polarized in such a 
way that a fraction £ is transverse while a fraction (1 — £) is longitudinal. The transversely 
polarized component of the B* gives rise to an angular distribution of (1 + cos 2 9) /A weighted 
by the relative probability £, while the longitudinal polarized component has an angular 
distribution of (1 — cos 2 9)/ 2 weighted by the relative probability (1 — £), where 9 is the angle 
between the outgoing photon 3-momentum and the polarization axis in the B* rest frame. 
We define the polarization axis to be the direction of the 3-momentum of the B* in the 
laboratory frame. Summing over all the helicity states of the B*, the angular distribution of 
the emitted photon is given by 

-, n ~ 1 + — t cos d , ( 4 ) 

dcos9 V 2 - W 



= 1 - 




(5) 



(6) 



If the relative probability T/L equals 2 as suggested by the naive spin counting rule, £ equals 
2/3 and the resulting decay angular distribution will be flat. In other words, there is no spin 
alignment of the B* meson. We will show that the above isotropic scenario is only true in 
the heavy quark mass limit and this limit is broken by the finite charm and bottom quark 
masses. 

This paper is organized as follows. In Section II, we will derive the polarized fragmenta- 
tion functions D^ B , (z, s) that depend on both the usual fragmentation variable z and the 
variable s which measures the virtuality of the fragmenting b quark. Covariant expressions 
for the transverse and longitudinal polarization sums for massive spin 1 objects will be de- 
rived so that covariant calculation can be performed. In Section III, we will discuss the heavy 
quark mass limit of mj ^> m c ^> Kqcd- We show that our polarized fragmentation functions 
satisfy heavy quark spin symmetry in this limit, and that heavy quark symmetry breaking 
arises from the next-to-leading term in the heavy quark mass expansion. In Section IV, we 
will study the spin asymmetry parameters £, a, and W in detail. We will point out that 
the anisotropic angular distribution in the two body decay of the excited meson (D*, B*, or 
B*) first arises at the next-to-leading term in the heavy quark mass expansion. We will also 
introduce the z-dependence in the spin asymmetry parameters £(z), a(z), and W(z), and 
study their variations with z. In Section V, we will discuss the scaling behaviors of the mean 
longitudinal momentum fraction (z) for the B* meson, for quarkonia, and for the heavy-light 
excited mesons. In Section VI, we will present the results of the fragmentation functions 
that depend on both the variable z and the transverse momentum pj_ of the B* meson with 
respect to the fragmentation axis, which is defined as the direction of the 3-momentum of 
the fragmenting b quark in the laboratory frame. Spin asymmetry parameters depending on 
p_L are also introduced and studied in detail. We conclude in Section VII. 
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II. Fragmentation Functions for Polarized B* Meson 

The derivation of the polarized fragmentation functions ofb—>B* follows closely to the 
unpolarized case 0, but requires separate contributions from the longitudinal and transverse 
components of the B* meson. For the unpolarized case all the helicity states are summed 
by the following formula 

E^(p)^) = -<r + ^, (t) 

A 

where p and M are the momentum and mass of the meson, respectively. Since we are working 
entirely within the nonrelativistic approximation for the heavy quark bound state, we will 
take M = + m c . The two transverse helicity states are usually summed by using the 
following equation 

E#(PK(p) = **-|jp (8) 
This noncovariant expression for summing the transverse polarizations often makes the ma- 
nipulation unnecessarily cumbersome. Here we present covariant formulas for summing the 
longitudinal and transverse helicity states of massive spin 1 objects. Recall that the longi- 
tudinal polarization 4- vector e^(p) is usually written explicitly as 

4« -(£•!&)■ « 9 > 

where E 2 = M 2 + \p\ 2 . Let us define an auxiliary 4- vector = (1, —p/\p\) such that n 2 = 
and n ■ p — E + With the help of this auxiliary vector one can rewrite e£(p) in the 
following covariant form, 

4.(p) = j7- — ■ (io) 

M n ■ p 

We can then obtain the following covariant expressions 

eTipKip) = KC~ — (p"*" + P"*") + l^-T^n" , (11) 
M z n ■ p [n ■ p) 

for the longitudinal polarization sum, and 

1 M 2 

^T(pVt(p) = ST + (p^ + pW) - -n»n\ (12) 

y n ■ p [n ■ p) 1 



for the transverse polarization sum. In general, these covariant formulas are not only useful 
in our calculation but also in many other applications. 

The initial fragmentation functions D^ B „(z, /x ) and Dj_^ B „ (z, /i ) are obtained by mod- 
ifying the calculation in Ref. || using Eqns. (|TT|) and (|T2| ) to separate the contributions 
from the longitudinal and transverse components of B*. The starting point is the following 
expression || 

r ( M 2 m 2 \ 
D^ B .{z) = JdsOls- — - D- b _+ B *(z, s) , (13) 

with 

1 \M\ 2 
D- b ^ m (z,s) = lim — — . (14) 

lD7r q /m b ->oo \JVlo\ z 

In Eqn.(|l4D, M. is the amplitude for producing a B* and a c quark from an off-shell b* with 
virtuality s = q 2 , where q is the 4- momentum of the b quark (the leading contribution is 
an one-gluon exchange diagram given in Fig. P; and .Mo is the amplitude for producing an 
on-shell b with the same 3-momentum q. If the momentum of the B* is p M = {po,Pi,P2,P3), 
in a frame where q^ = (90,0,0,93), the longitudinal momentum fraction z is defined by 
z — {.Po + P3)/{Qo + ^3)- The kinematical relation among the three variables z, s, and 
p± = \p±\, where p± = (pi,P2) is the transverse momentum of the B* meson with respect to 
the fragmentation axis, is 

_ M 2 + p\ | m 2 c +p\ 
z 1 — z 



Note that the 6 function constraint in Eqn. ([13| ) arises from the positivity of p±. 

The amplitude A4 is evaluated in the axial gauge with an auxiliary 4- vector defined 
above in the covariant polarization sums. In this gauge, the dominant contribution arises 
from the Feynman diagram depicted in Fig. [I]. Other diagrams are suppressed by powers 
of mfc iC /g - In other words, factorization is manifest in this gauge 0. The amplitude M. is 
given by 



/ iVc m c (s — ml) 2 



.s — ml 



xt I 2M{i + m b )f x {p) + n .( q _ rp) W + M)ft(p)i\ v(p') , (16) 



where p and p' are the momentum 4- vector of the outgoing B* and c respectively, e\(p) is 
the polarization 4-vector of the B* meson with helicity A, -R(O) is the radial wavefunction of 
the (be) bound-state at the origin, M = m& + m c is the mass of the bound state, r = m c /M, 
f = 1 — r, and T is a symbolic Dirac structure representing the source to create the energetic 
b quark in the hard subprocess. In Eqn. (0), Cp = (N 2 — 1)/2N C , where N c is the number 
of color. -R(O) can be determined by a potential model calculation |T1 or extracted from the 



B c decay constant /b c . The latter can be calculated on a lattice or measured in the future 
experiments. The relation between -R(O) and fs c is given by |i?(0)| 2 = 7rMs c /^ c /3. Squaring 



the amplitude A4, summing over the color and spin of the c quark, and using Eqn. (11) to 
project the longitudinal component of the B* meson, we get 

M 5 1 

£ \M\ 2 = 7ia 2 s C F \R(0)\ 2 — 7 2^Tr(rr^) x A L (z, s) , (17) 

m c [s — m b ) 

where A L (z, s) is given by 

A L (z,s) = —(l-2fz + f(l-2r)z 2 ) 

+ ~ m ',L f-4 + 2(3 - Ar)z - (1 - 8r + Ar 2 )z 2 - f(l - 2r)z 3 ) 



+ 



z(l - fz)M 2 
4(1 - z)(l + rz) 2 (s - ml) 2 



(1 - fz) 2 M 4 ' y ' 

To arrive at this expression we have made the substitutions of p = zq and p' — (1 — z)q at 
the final step, which are accurate to leading order in mb, c /Qo- I n the fragmentation limit of 
qo/rrib — > oo, the tree level amplitude A4o is simply 

^|^o| 2 = iV c Tr(rr^). (19) 



Thus we obtain 



where we have defined 



_ a 2 s (2m c )C 2 \R(0)\' - );| 
2 l.V, ~///; ! 
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The scale of the strong coupling constant has been set to be 2m c - the minimal virtuality of 
the exchange gluon ||. Doing the s integration, we obtain the longitudinal fragmentation 
function 



D 



0,/i )= n 



rz(l — z) 



2 + 2(2r - 3)z + (16r 2 - lOr + 9)z 2 



(1 - fzf 

- 2f(6r 2 -5r + 4)2 3 + 3f 2 (2r 2 -2r + l)^ 4 



(22) 



where the initial scale \x§ has been set to be (m^ + 2m c ) - the minimal virtuality of the 
fragmenting b quark ||. 

Similarly, one can use Eqn. (0) to project out the transversely polarized B* state and 
deduce 

A T (z,s) 



Di^ B *{z, s) = -rNM e 



(s - m 2 b ) 



2N4 ' 



(23) 



with 



A T (z,s) 



16 , 

— 1 - 2fz + f(l +r)z' 



8(s — mT) , , s 9 o 

+ — t~ — _ x \' rn (2 + 2(r - 2)z + (2r + l)z 2 + fz 3 



+ 



*(1 - fz)M 2 
8z 2 (l- z){s-m 2 b ) 2 
(1 - fz) 2 M 4 



(24) 



Doing the s integration, we obtain the transverse fragmentation function 

rz(l — z) 2 



D 



b^B* 



2N 



'1-fzY 



2 + 2(2r - 3)z + (10r 2 - 4r + 9)z 2 
-2f(r + A)z 3 + 3fV 



(25) 



At this point, a few cross checks can be made. Adding -Df^o* (z) and Dl_^ R * (z), we reproduce 



the unpolarized result given in Refs. ||, |10|- By setting r = 1/2 in the Eqns. ( f22j) and (25), we 
reproduce the results of c — ► polarized given in Ref. [0. The results given in Eqns. ( f22|) 



and pq) also agree with a recent calculation of Ref. [15 



One can easily extend these results to the case where the initial fragmenting b quark is 
also polarized. Let us denote D^\)(z) to be the fragmentation function for a heavy quark 
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Q with helicity h = ±1/2 to split into a vector meson V* with helicity A = 0, ±. Thus, by 
definition, we have 

D L {z) = \(D {hQ) {z)+D { _ hQ) {z)) , (26) 

and 

D T {Z) = i (£>(!.+)(*) + + + • (27) 

Parity invariance implies 

D {h;X) (z) = D^ x) (z) . (28) 
Therefore, we deduce the following relations 

D ih0) (z) = D i _ h0) (z)=D L (z) , (29) 

and 

+ = D { _ h+) {z) + U ( _ iH («) = D T {z) . (30) 

These relations immediately imply that the polarized heavy quark fragmentation functions 
D L (z) and D T (z) are the same whether the initial heavy quark is polarized or unpolarized. 



III. Heavy Quark Spin Symmetry 



Using the technique of HQET |T], 0, [| , Jaffe and Randall |16j have recently shown that 



the fragmentation functions Dq-^h(z) for a heavy quark Q to split into a hadron H with 
one heavy constituent quark can be expanded as a power series in r, 

D Q ^ H (z) = ^- + b(y) + 0(r), (31) 

where a(y) and b(y) are functions of y = (1 — fz)/rz and C(r) denotes all other terms 
higher order in r. The leading term a(y) is independent of the heavy quark spin and flavor; 
while the next-to-leading term b(y) and all higher order terms contain heavy quark spin- 
flavor symmetry breaking effects. One can verify easily that our polarized B* fragmentation 
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functions can be expressed in this form by carefully expanding the powers of r and (1 — rz). 
The results are 



N(y 



Dl 



y b 

2N(y 



1 



+ 4y + 3y 2 )~(8-8y + 5y 2 + y 3 ) + 



If 



-(8 + Ay + 3y 2 ) - (8 - 8y - y 2 + y 3 ) + .. 



(32) 
(33) 



Obviously, the leading terms of D L (z) and D T (z) in Eqns. ( |52"D and ( |3"3"D differ only by a factor 
of 2 and thus obey heavy quark spin symmetry. The O(r ) terms and beyond in Eqns. ( |52"D 
and ( |3"3"D are different due to heavy quark spin- flavor symmetry breaking effects. In fact, one 
can show that |T7| the leading order terms in Eqns. ( |3"2"| ) and fl3"3] ) can be derived by using 
the Feynman rules of the leading operator in the HQET Lagrangian, while the O(r ) pieces 
arise not only from the next-to-leading (1/M) operators but also from the small component 
of the heavy quark spinor in the HQET. 



In Fig. |2|, we plot the full PQCD fragmentation functions D 



z) and DT D *(<2) from 



'b-*B, 



Eqns. ( p22[) and (p25|) , and the corresponding heavy quark mass expansion from Eqns. ( p2|) 
and (|33l) . We take mt, = 4.9 GeV and m c = 1.5 GeV throughout the paper. We employ 
a simple form of a s by evolving from its well measured experimental value at the Z°-mass, 
namely 



a s (n) 



a s {m z ) 



(34) 



1 + (b/2n)a s (m z )\og([i/m z ) ' 
where b = (lliV c — 2n/)/3, rif is the number of active flavors at the scale n, and a s (mz) = 
0.12. It is clear from Fig. [| that the sum of the leading and the next-to-leading terms in the 
heavy quark mass expansion is a very good approximation to the full PQCD result, and the 
difference is of order 0(r). We note that the width between the peak and the endpoint of 
the fragmentation functions scales as r. 



Since our polarized fragmentation functions given in Eqns. (22) and (25) are consistent 



with the general analysis of Jaffe and Randall [16], we can apply them to describe the 



heavy-light mesons as well. The formulas given in Eqns. (^) and (25) can be regarded as 
phenomenological fragmentation functions with two free parameters iV and r, to describe 
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the nonperturbative process of a heavy quark splitting into a polarized heavy-light excited 
meson, such as c — > polarized D* and b — > polarized B*. The parameter N can be adjusted 
to describe the overall normalization and r is the mass ratio of the light constituent quark 
mass to the meson mass. Although in the D* and B* systems there are probably large 
nonperturbative and relativistic effects that we have not taken into account, our perturbative 
QCD fragmentation functions with the free parameters N and r can at least provide some 
insights to these systems while precise nonperturbative fragmentation functions for c — > D* 
and b — > B* are not available yet. Our PQCD fragmentation functions only depend on two 
free parameters as does the phenomenological Peterson fragmentation function [ljj, which 
is widely used in the literature to describe c — > D* and b — > B* . However, our fragmentation 
functions also carry spin informations. We expect our functions should be equally successful 
as phenomenological description of fragmentation, but they are more predictive since they 
also contain spin informations. They also have an additional virtue of being rigorously 
correct in some limit, namely, the limit in which the lighter quark mass is much greater than 
Aqcd- 

Finally, we note that the 2-to-l spin counting ratio for the transversely and longitudinally 
polarized states is only true for the S-wave excited states in the heavy quark limit. For the 
P-wave excited states, this statement is incorrect M. 



IV. Spin Alignment in Heavy Quark Fragmentation 

In this Section, we will study the spin alignment of the excited S-wave heavy meson 
(J/ip, T, B*, D*, or B*) produced by heavy quark fragmentation. To begin with, we 
integrate D L ' T (z, /io) over z to get the total fragmentation probabilities, or equivalently the 
first moments D L,T {1) of the corresponding fragmentation functions: 



D L {1) = ( l dzD L {z^ ) 
Jo 



N 



24 + 89r - 486r 2 + 354r 3 + 289r 4 
15 f 5 
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r(7 - 16r - 9r 2 + 30r 3 + 6r 4 ) , . . 
+ ~ R log(r) 



D T {1) = f 1 dzD T (z,fx ) 
Jo 

"24 + 119r + 54r 2 + 84r 3 - llr 4 r(7 + 2r + 9r 2 ) 



2N 



log(r) 



(35) 



(36) 



where r is defined as the mass ratio mnght quark /"^meson, just the same way as we defined for 
the B* meson. From now on we will omit the subscript b — > B* in the fragmentation function 
and understood that it can refer to either B*, J/ip, T, or the heavy-light mesons D* and 
B* . It is well known that the first moment of fragmentation function has zero anomalous 
dimension and hence it does not evolve with the scale. We therefore drop the \x dependence 
in the first moments D L,T (1). 

Since the dominant production mechanism of S-wave excited meson states at the large 
transverse momentum region is due to fragmentation, we can identify the quantities L and 
T, defined in Section I, to be the first moments D L (1) and D T (1), respectively. The various 
spin asymmetry parameters introduced in Eqns. (|l|) - @ can be expressed as ratios of the 
first moments of the fragmentation functions as follows: 

e = DT{1) (37) 

? D T (1) + D L (l) ' K ' 

2D\\) -D T (1) 

a ~ — mi) — ' ( } 

W ~ ZF(1)+2^(1) • [6J) 
Note that these spin asymmetry parameters depend only on the parameter r and do not 
depend upon the overall constant N and the evolution scale \i. 

Before we proceed further, it is instructive to repeat here the usual arguments of how the 
heavy quark spin information is lost during hadronization of the heavy quark into a S-wave 
heavy-light meson. Suppose a spin down heavy quark combines with a spin up or spin 
down light anti-quark q forming the state Q^(f or Q^-qK Since parity is conserved in the 
fragmentation process, these two states must occur with equal probability. While Q^g 1 is an 
eigenstate of the total spin S = 1; Q^(f is a mixture of the spin states S = and 5 = 1. 
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One can then decompose the state Q^-q^ into a sum of eigenstates of the total spin S = 
and S = 1: 



^ = 75 



1 

+ 75 



v/2 



v/2 

These 5 = and S = 1 components are identified as the pseudoscalar P and vector meson 
V*, respectively. In the limit of rriQ — ► oo, these two states P and ^* are degenerate. They 
will then have the same time evolution and will propagate coherently. The spin wavefunction 
will remain the same as Q^q^ and the coherent superposition of the two meson states will 
preserve the heavy quark spin forever. In reality, mq ^ oo; the pseudoscalar P and vector 
meson V* have a slight mass difference AM and hence different time evolutions. In most 
cases, like the heavy-light mesons and the (be) mesons, the finite mass difference AM is 
considerably larger than the decay rate Ty* of V* — > P + X, where X denotes a photon 
or a pion. At a time t ~ AM -1 , the S = and S — 1 components become completely 
out of phase and incoherent before any decay actually occurs. Thus the heavy quark spin 
is depolarized over a period of time characterized by the chromo-magnetic dipole moment, 
which is also responsible to the finite mass difference AM. At a later time t ~ T v l, the 
vector meson V* decays into the pseudoscalar P. The spin information of the heavy quark is 
carried in the relative population of the helicity +1 and helicity —1 components of the vector 
meson. If the decay proceeds through electromagnetic or strong interactions, the parity 
invariance combined with rotational symmetry implies that the spin information is lost in 
the angular distribution of the final state particles. Nevertheless, the angular distribution 
carries important information about the fragmentation process through the spin alignment 
variables, which measure the relative population of the helicity and helicity ±1 components 
of the vector meson. 

In the heavy quark mass limit of mj, 3> m c ^> Aqcd, t — > 0, the spin asymmetry 
parameters obtained in Eqns. (|37D -(|39D take the following simple forms 

e-| = ^ + 0(r 2 ), (41) 
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a = ~r + 0(r 2 ), (42) 

2W - 1 = lr + 0(r 2 ) . (43) 
8 

Note that, to order r, the subleading terms 0{r\ogr) in Eqns. fl3"5]) and ( |3"6D cancel in these 
quantities. These limits make it clear that the angular distribution given in Eqns. (f§)-(§|) 
is isotropic to leading order in r. This result is in accord with the recent general analysis of 
Falk and Peskin [^] using heavy quark symmetry. As a matter of fact, in the heavy quark 
limit the spin asymmetry parameter W coincides with the Falk-Peskin variable W\/2, which 
is the conditional probability for a heavy quark fragmenting into a heavy-light meson with 
a spin | light degree of freedom to be in the helicity states hi = +| or — ~. Light degree of 
freedom denotes collectively the light quark plus all the soft gluons that combine with the 
static heavy quark to form the color singlet heavy-light system. Parity invariance implies 
these two helicity configurations (hi = ±|) of the spin \ light degree of freedom must occur 
with equal probability J7J. Thus 2W — > Wi/ 2 = 1 in the heavy quark limit. Consequently, the 
anisotropy of the decay product of the excited S-wave heavy-light meson arises entirely from 
the spin-flavor symmetry breaking effects of the heavy quark. Our results of the polarized 
fragmentation functions and the spin asymmetry parameters obtained above allow us to 
study these symmetry breaking effects as a function of r. 

We plot these three spin asymmetry parameters versus r in Fig. ^|(a), (b), and (c) re- 
spectively. The anisotropy of 5.7% (T^j in the decay distribution of J/tj) — > ££, where £ = e~ 
or fi~, can be immediately evaluated from the values of £, a, or W at r = 1/2 in Figs. [3[ 
It means that the decay distribution behaves like (1 + 0.057 cos 2 6), where 9 is the angle 
between the outgoing lepton and the polarization axis of the J/ip. In principle, a 5.7% 
asymmetry is measurable, but it is seriously contaminated by the J/ip's coming from B 
decays. This anisotropic distribution is also true for T — > It. Surprisingly, the anisotropy in 
B* — ► Bc'-f is only ~ 5.8%, almost the same as in the equal mass quarkonium case. But it is 
relatively clean because the production rates for the 2S, 2P, 3P, and 3D (be) states which 
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can contaminate the direct b — > B* fragmentation by their hadronic cascades or radiative 
decays into B* are in general small. The production rates for the 2S (be) states are about 
60 % of the corresponding IS" states l|. The fragmentation probabilities for b — > P-wave 



(be) states are only about 10 % of the S-wave case [19|. While the fragmentation functions 
for b — > D-wave (be) states are not known, they are not expected to be large. In order to 
measure the asymmetry, one still has to disentangle the issue of whether the B* is coming 
via direct b fragmentation or cascade from higher excitations. The anisotropies in the D* 
and B* systems are 5.1% and 2.6% respectively. 

Since the ^-integrated spin asymmetry parameters in general imply small anisotropies, 
one might try to get more sensitivity by studying the ^-dependent spin asymmetry parame- 
ters defined by: 

These ^-dependent spin asymmetry parameters are necessarily \i dependent because the 
fragmentation functions depend on the scale \i. They are shown respectively in Fig. f|(a), 
(b), and (c) at the corresponding initial scale \iq and r for J/tp(r = 0.5), B*(r = 0.23), 
D*(r = 0.17), and B*(r = 0.058). We have chosen /io as the sum of the heavy constituent 
quark mass and double the light constituent quark mass for each bound state, and r to be 
the ratio of the light constituent quark mass to the meson mass. We also take the light u or 
d constituent quark masses in the B* and D* mesons to be 0.3 GeV. The dependence of £(z), 
ct(z), and W(z) on z in Fig. [| shows that the maximum asymmetry occurs at z around 0.6 - 
0.8. We note that despite having a different initial scale /io, the spin asymmetry parameters 
for the T are the same as those for the J/if>. 

One might worry about how the scale \i affects the shapes of £(z, fi), a(z, /i), and W(z, fi). 
Here we remind the readers that the z-integrated £, a, and W are independent of the scale. 
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We can evolve the polarized fragmentation functions by solving the Altarelli-Parisi evolution 



equation [20| 

^D L ' T (z, p) = f 1 P„(£) ^(y, a*) , (47) 

^ J* y v 

where P qq {x) = ^^-Cf[(1 + x 2 ) /(l — x)] + is the usual quark-quark Altarelli-Parisi splitting 
function, and the plus distribution is defined by f(x) + = f(x)—5(l—x) Jq dx' f(x'). When the 
initial heavy quark is polarized, we should use the polarized Altarelli-Parisi splitting function 
AP qq (x), which happens to be the same as P qq (x), in Eqn. (f|7|). However, in the small z 
region, the Altarelli-Parisi evolution equation does not handle the threshold effect properly 
f2lfl . Unphysical behaviors occur at the small z region for the fragmentation functions and 
hence for the spin asymmetry parameters, as one evolves the scale up to, say, fi = mz/2- 
Therefore, we can only discuss the evolution behaviors at the large z region with confidence. 
We found that in the large z region the curves for the spin asymmetry parameters at the 
scale fi = mz/2 do not differ significantly from those given in Fig. |4], both in shape and 
magnitude. In the future when detailed fragmentation data for the polarized D* and B* 
become available at LEP, it will be very interesting to compare our theoretical predictions 
of the spin asymmetry parameters in the large z region given in Fig. |] to the experimental 
results. 

A recent set of experimental data on the spin alignment of D* meson was from the 
CLEO detector operating at a center-of-mass energy of 10.5 GeV p2| . The data set given 



in Table I of Ref. [22] is for the spin asymmetry parameter a(z). The scale of the data 
set is taken to be half of the center-of-mass energy, i.e., 5.25 GeV, which is not very far 
from the /j,q = 2.1 GeV of the D* meson. We therefore ignore the evolution effects and 
directly compare the D* curve in Fig. £|(b) with the data from the CLEO measurements. 
Unfortunately, since the experimental data for a(z) had very large error bars, we have to 
show the comparison on another graph (Fig. [5D with a larger vertical scale. 

While our model always predicts a slightly negative value of a(z) for all values of z, 
the CLEO data are rather scattered with statistical tendency toward the positive side, as 
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indicated by the mean value (oc(z)) = 0.08 ± 0.07(stat) ± 0.04(sys). As shown in Fig. |5], the 
agreement is good because our curve a(z) is within 2a of all the data points. In the CLEO 
analysis, they also concluded that the data was marginally consistent with zero. 

In this Section, we have shown that both the ^-integrated and the ^-dependent spin 
asymmetry parameters indicate a small misalignment in the two polarization states of the 
excited S-wave heavy mesons produced by heavy quark fragmentation. As indicated by the 
parameter a(z), one predicts that the transverse states produced by heavy quark fragmen- 
tation should be populated slightly more than would be given by naive spin counting over 
the entire physical region of z. 



V. Mean Longitudinal Momentum Fraction 



Another useful observable, the mean longitudinal momentum fraction (z), of the B* (or 
the heavy- light mesons D* and B*, or the heavy quarkonia) at the scale /i is defined as 

Jo 1 dzzD(z, fi) 



Jo dzD(z,n) 



(4* 



In other words, (z) is the ratio of the second to the first moment of the fragmentation function 
at the scale /i. Since the anomalous dimensions of all the moments of the fragmentation 
function are known explicitly, the scaling behavior of (z) can be determined as 



£(2,/i) 
D(l) 



<*>c 



Q!s(A*o) 



-■ -22 



(49) 



where (z) = D(2, fx )/D(l), 7 = -4C F /3, and b = (11N C - 2n/)/3. The first moments 
D L,T (1) of the longitudinally and transversely polarized fragmentation functions are given 
in Eqns. (j3~5[) - (j36|) . The second moments D L ' T (2, /i ) at the scale /j,q are given by 



D L (2,v ) = f 1 dzzD L (z,fi ) 
Jo 





2N 



12 + 97r - 453r 2 + 252r 3 + 437r 4 + 15r £ 



+ 



15f 6 

r(5 - llr - 15r 2 + 33r 3 + 12r 4 ) 



log(r) 



(50) 
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D T (2,fi ) = [ 
Jo 



1 



AN 



dzzD (z, /i ) 
12 + 112r - 33r 2 + 252r 3 + 17r 4 



r(5 + r + 9r 2 + 9r 3 ) 



log(r) 



15f 6 



+ 



(51) 



In Fig. |, we plot (z) L ' T versus the scale \i for the four meson systems J/ip(r = 0.5), B*(r = 
0.23), D*(r = 0.17), and B*(r = 0.058) with the same input parameters defined in the pre- 
vious Section. Due to slightly different value of a s (/i ), the corresponding curves of (z) L ' T 
for the T(r = 0.5) are slightly different from those of the J/ip(r = 0.5). We will not present 
the T curves here. 

As the scale fi increases across each heavy quark threshold (2m c and 2m;,), the number rif 
of active flavors increases by one unit. The kinks on the curves at /i = 2m c and Im^ in Fig. |6] 
are due to these threshold effects. Notice that only slightly noticeable differences between 
the longitudinal (z) L and transverse (z) T occur for the J/ip and B*, but no differences can 
be seen for the D* and B* mesons. Hence, experimentally using (z) L,T to distinguish the 
polarizability in heavy quark fragmentation into the S-wave excited meson states is not 
plausible. We list the {z)q at \i = //o and (z) at \i = rriz/2 for the four different mesons 
in Table [l|. Numerically, there are no noticeable difference between the longitudinal and 
transverse polarizations for all the four meson systems. Therefore, to a good approximation, 
one can set these polarized (z) L,T values for each meson to be the unpolarized (z) value of the 
corresponding meson. Experimentally, unpolarized (z) are available from the LEP, CLEO, 
and ARGUS data. The measured quantities are (xe) c -^d*, ( x e)c->h c , and (xE)b->H b , where 
xe is the energy fraction of the hadron or meson relative to one half of the center-of-mass 
energy of the machines, and H c and denote the charm and bottom hadrons, respectively. 
Because xe is a good approximation to the fragmentation variable z, we simply treat them to 
be the same in the following comparisons. Since the data for c — > D* is available separately, 
we will use (xe) c ^d* instead of the inclusive (xe) c ^h c - O n t ne other hand, only the inclusive 
(xE)b^H b has been reported at the LEP. Nevertheless, this inclusive value of {xE)b^H h should 
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be close to the {xE)b^B*, since b — > B* + X is expected to be the dominant fragmentation 
mode of the 6-quark. In addition, we take the scale of the various measurements to be one 
half of the center-of-mass energies of the machines. 

Next we will describe briefly how we obtain the average from the LEP, CLEO, and 
ARGUS data. For the LEP data measured values of {xe) c ^d* are from OPAL (0.52±0.0316), 
ALEPH (0.504 ± 0.0188), and DELPHI (0.487 ± 0.0158) g|, in which we have already 
combined their systematic and statistical errors in quadrature if they are given separately. 
We then simply take the mean of the central values from the three experiments to be the 
average central value. For the combined error we add the absolute errors from the three 
experiments in quadrature and divide it by 3. We thus obtain (xe)c^d* = 0.504 ±0.0133 by 
combining all three LEP data. Similarly, we have measured values of (xE)b^H b from OPAL 
(0.726 ± 0.023), ALEPH (0.67 ± 0.050), DELPHI (0.695 ± 0.0326), and L3 (0.686 ± 0.017) 
p4| . Repeating the same exercise, we obtain {xE)b->H b = 0.694 ± 0.0166. 

The data from CLEO [^5] and ARGUS [2B| were given in form of fragmentation functions. 
We need to calculate (xe)c^d* from their fragmentation data. In the CLEO paper, the 
fragmentation function of c — > D* were given in Table I (a) and (b) of Ref. JHJ, which 
correspond to two different detection channels. We combine the two tables with the y-value 
of each bin equal to the mean of the two, and the error of the y- value of each bin equal to one 
half of the two errors added in quadrature. Then the value (xe) is obtained by evaluating the 
two integrals in the ratio Jq 37 cIxeXeD(xe) / J0.37 dxED(xE) by the method of discrete sums. 
Assuming the error is only in the y- value of each bin, we obtain (xe) c -^d* — 0.654 ± 0.0563, 
where the error is obtained by adding the errors from each bin in quadrature. Similarly, 
from the ARGUS data we obtain (xe) c ~^d* — 0.642 ± 0.067. Finally, we combine the two 
averages from CLEO and ARGUS and obtain (x E ) c ^ D * = 0.648 ± 0.0438. Since the CLEO 
and ARGUS operating center-of-mass energies were so close to each other (10.55 GeV for 
CLEO and 10.6 GeV for ARGUS), we assume the scale of both measurements to be 5.3 GeV. 

The LEP average for {xe) c ^d* and {xE)b~*u b at /i = m^/2, as well as the combined 
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CLEO and ARGUS average for (xe}c->d« at /i = 5.3 GeV are shown in Fig. ||. Excellent 
agreement between our predictions and the data is demonstrated. Here we remind the readers 
that for the B* and D* mesons, we have assumed a nonrelativistic bound state picture with 
the light constituent quark masses being set to be 0.3 GeV. Uncertainties arising from the 
overall normalization of the fragmentation functions are cancelled, as indicated by Eqn. (pE8|). 
Therefore, the uncertainties in the heavy-light meson systems come only from how we define 
the strong coupling constant a s , the initial scale Hq for the fragmentation functions, and what 
values we choose for the light constituent quark masses. Summing up, the mean longitudinal 
momentum fractions {z) L,T for a heavy quark to fragment into polarized vector mesons do 
not show any measurable difference between the longitudinal and transverse components. 
Nevertheless, the predictions of (z) by our fragmentation functions for D* and B* mesons 
at different scales are in excellent agreement with the measured data from the LEP, CLEO, 
and ARGUS detectors. 



VI. Transverse Momentum Dependence 

Sections III and IV showed that the ratio of longitudinal-to-transverse populations from 
6-quark fragmentation into polarized B* mesons is only marginally different from that given 
by the naive spin counting (i.e. the heavy quark mass limit), as indicated by both the 
^-integrated and z-dependent spin asymmetry parameters (see Figs. |3] and |]). Similar con- 
clusions hold for the other mesons J/ip, T, D*, and B*. As a consequence, the two body 
decays of the excited S-wave mesons into the corresponding pseudoscalar ground states with 
emission of photon or pion are almost isotropic. In addition, we showed in Section V that 
it is not feasible to use the mean longitudinal momentum fractions (z) L ' T to distinguish the 
longitudinal and transverse polarizations of the excited S-wave mesons produced by frag- 
mentation of heavy quarks (see Fig. ^j. Hence, information about the spin of a heavy quark 
is not easy to extract from the fragmentation data of the heavy quark into S-wave excited 



mesons 27 . 
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So far, we have only investigated the z-dependence of the fragmentation functions and 
of the spin asymmetry parameters. All the dependence on the motion of the meson perpen- 
dicular to the fragmentation axis has been integrated out. In this Section we will investigate 
the dependence of the fragmentation functions on the p± = \p±\, where p± = (pi,P2) is the 
transverse momentum of the meson with respect to the fragmentation axis. Recall that from 
Eqn. ([15]) , we have 

p\ = z{\ - z )[s- — - ^-J . (52) 

Introducing the dimensionless variable t = p±/M and trading the variable s to t from 
Eqn. (|52|) , we can define the fragmentation functions D(z,t) and D(t) according to 



/*oo rl roc 

/ dtDit) = dz dtD(z,t) 
Jo Jo Jo 



dz ds 8 ( s 



z 



m. 



D(z,s) 



(53) 



This implies 



D{z,t) 



2MH 
z(l-z) 



-D(z,s) , with s = M 



1+t 2 



+ 



r 2 + t 2 
1-z 



(54) 



To our knowledge, QCD evolution equation for the fragmentation function D(z,t) that de- 
pends on both the longitudinal momentum fraction z and the rescaled transverse momentum 



t has not been written down. But formalism, like those in Refs. [pg], which dealt with similar 



issue in the parton distribution functions, may apply for fragmentation functions as well. For 
the following we will ignore the issue of QCD evolution effects in the fragmentation functions 
with explicit t-dependence. 



Integrating over z in Eqn. (54), we obtain the t-dependent polarized fragmentation func- 
tions D L ' T (t): 

Nr 1 



2f 6 t 6 (l + t 2 ) 5 



24fV(l + r) 2 log(r 



+ 12t{l + t 



2\2 



Ar 



r(2 + r + 2r 2 )t 2 + f 2 t 6 
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2\2 



+ 3(1 + t 2 ) 
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arctan 



rt 
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(55) 
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(56) 



As a cross check, one can integrate D L ' T (t,fi ) over t from to oo, and get back the total 
fragmentation probabilities given in Eqns. (55|) and (|36|). 

A couple of asymptotic behaviors of D L ' T (t, /i ) are in order. As t — > 0, D L ' T (t, Ho) 
vanishes linearly in t: 



D L (t,v ) 
D T (t,fi ) 



2Nrt 



6 log(r) + -(8 - 46r + 160r 2 + 101r 3 - 13r 4 ) 

35H 
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(57) 
(58) 



As t — > oo, D L,T (t, /io) falls off like 1/t 3 according to 
Ar 



D L (t,fJ, Q ) 



£ T (t,/i ) 



f 6 t 3 
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f 6 t 3 



12r(2 + r + 2r 2 ) log(r) + f(l + r)(6 + 20r + r 2 + 4r 3 - r 4 ) 
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The dependence of the fragmentation functions Dj^ B * it, /io) on t at the initial scale /io 
are shown in Fig. |7]. Both functions peak at t ~ 0.2, i.e., ~ 1.3 GeV for B* meson. 



(59) 
(60) 



24 



At large p±, one sees that the transverse D T (t) falls off more rapidly than the longitudinal 
D L (t). Explicitly, the curves for the longitudinal and transverse polarizations cross over at 
t ~ 1.5, i.e., p± ~ 10 GeV for B* mesons. This is a very clean sign of difference between 
the longitudinally and transversely polarized states of B*. Overall, the transverse states 
are populated about twice as much as the longitudinal one, but the longitudinal component 
becomes dominant beyond the cross-over. Unfortunately, it might be hard to observe the 
cross-over experimentally because the cross-over lies well beyond their peak values. Both 
the longitudinal and transverse spectra have dropped substantially by the time they reach 
the cross-over, which implies very small cross sections for p± > 10 GeV. In addition, there 
are large uncertainties in determining the fragmentation axis of a jet and consequently also 
the value of p± for the B*. Nevertheless, in principle, imposing a high p± cut can eliminate 
a large sample of the transverse population of the B*. 

We can also define the average transverse momentum (p±) by 

(p±) _ M _ Io Q °dttD(t) 

M W IS°dtD(t) ' 1 ' 

Numerically, the average (t) for the longitudinal and transverse B* are 1.1 and 0.61 respec- 
tively, and so the corresponding average (p±) are 7.0 GeV and 3.9 GeV, respectively. Unlike 
(z), there is a big difference in (pj_) between the longitudinal and transverse states of B*. 
The cross-over, and the average (t) and (p±) for the B*, J/ip, D*, and B* are all summarized 
in Table [2| Entries for T are the same as those for J/ip, except that the average pj_ for T 
is larger by a factor of mj,/m c . A significant difference in (t) between the longitudinal and 
transverse states persist in J/ip, T, D*, and B* systems. 

Despite the difficulties in their experimental measurements, it is instructive to evaluate 
the t-dependent spin asymmetry parameters defined by 

«" = wwrkf)' (62) 

2D L (t) - D T (t) , . 

«(*) = )j T{t) , (63) 
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In Fig. we plot the three spin asymmetry parameters a(t), and W(t) as functions of t 
for the mesons that we are considering. The curves in Fig. [5] vary far more dramatically than 
the corresponding curves in Fig. |], because the longitudinal and transverse fragmentation 
functions show a very different dependence on p±. For all the vector meson systems that we 
are considering, the t-dependent longitudinal and transverse fragmentation functions cross 
over at t ~ 1.5 — 2.0. Consequently, the corresponding spin asymmetry parameters also 
change very rapidly at t ~ 1.5 — 2.0. 

With D L ' T (t) replaced by D L,T (z,t) in Eqns. (|62"D - (|64|) , one can also introduce the spin 
asymmetry parameters £(z,p±), a(z,p±), and W(z,p±) that depend on both the variables z 
and p±. It may be possible to measure such spin asymmetry parameters, we will not discuss 
them any further here. 

VII. Conclusions 

In this paper, we have used PQCD to derive the polarized fragmentation functions to 
leading order in a s for a b antiquark fragmenting into the longitudinally and transversely 
polarized B* states. Parity invariance also implies that the polarized heavy quark fragmen- 
tation functions D L (z) and D T (z) derived in this paper are the same whether the initial 
heavy quark is polarized or unpolarized. These polarized fragmentation functions can be 
used to define various spin asymmetry parameters that can be determined experimentally. 

In this work, we have also used the PQCD fragmentation functions to study the spin align- 
ment of the other S-wave excited mesons that carry charm or beauty. The spin asymmetry 
parameter a(z) for c — > D* is consistent with the measurements by the CLEO detector. 
The mean longitudinal momentum fractions (z) predicted for D* and B* are also in ex- 
cellent agreement with the measurements of the LEP, CLEO, and ARGUS detectors. We 
also demonstrate a very interesting dependence of the polarized heavy quark fragmentation 
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functions and the spin asymmetry parameters on the transverse momentum p± of the vector 
meson relative to the fragmentation axis. Longitudinally polarized vector mesons have a 
harder p± spectrum than the transversely polarized states. 

The perturbative QCD fragmentation functions that we obtained in this paper and in 
Refs. |^9], H are expected to work well for c — > r) c , J/ip, b — > r] b , Y, and b — > B c , B* in any 
high energy processes with large transverse momentum px- They are rigorously correct in 
the limit that the heavy quark masses m c and mj are much larger than Aqcd- Corrections 
to the fragmentation contributions are of order M 2 /p T (M is the mass of the meson) and 
therefore small at large enough px- 

Our fragmentation functions also seem to provide a successful phenomenological model for 
describing the spin dependence of charm and bottom fragmentation into heavy-light mesons 
like D* and B* . Furthermore, our fragmentation functions are consistent with heavy quark 
symmetry. Dominant error comes from the neglect of relativistic corrections and higher order 
perturbative corrections. We conclude that the PQCD-inspired fragmentation functions 
could be useful in describing charm and bottom fragmentation in the future experiments. 
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Figure Captions 

1. The leading order Feynman diagram contributing to the fragmentation process b — > B*. 

2. The polarized fragmentation functions D^ B ,(z) and Dj_^ B ,(z) versus z at the initial 
scale fiQ. The sum of the first two terms in the heavy quark mass expansion are also 
shown for the longitudinal and transverse polarizations. 

3. The spin asymmetry parameters (a) £, (b) a, and (c) W versus the mass ratio parameter 
r. 

4. The spin asymmetry parameters (a) £(z,fj), (b) a(z,/i), and (c) W(z,fi) versus z at 
the initial scale //q for r = 0.5 (solid), 0.23(dash), 0.17 (dotdash), and 0.058 (dot) 
corresponding to J/ip, B*, D*, and B*, respectively. 

5. A comparison between the spin asymmetry parameter a(z) for c — > D* predicted by 
our model and the experimental measurements by the CLEO Collaboration p2] . 



6. The mean longitudinal momentum fraction (z) L,T versus the scale \i at r = 0.5, 0.23, 0.17, 
and 0.058 for the J/ip, B*, D*, and B*, respectively. The longitudinal curve is solid 
and the transverse one is dotted. The two measurements (xe) c ^d* and {xE)b->H b from 
the LEP detectors are shown at fi = mz/2, and the combined CLEO and ARGUS 
measurement (xe)c^d* is at [A — 5.3 GeV. 

7. D^ B „(t) and D^ Bt (t) versus the rescaled transverse momentum t = p±/M of the B* 
meson at the initial scale /j,q. 

8. The spin asymmetry parameters (a) (b) a(t), and (c) VV(t) versus t = p±/M at 
r=0.5 (solid), 0.23 (dash), 0.17 (dotdash), and 0.058 (dot), which correspond to J/ip, 
B* D*, and B* mesons, respectively. 
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Table 1: The mean longitudinal momentum fraction (z) when a heavy quark fragments into a 
polarized heavy meson in the spin-orbital state 3 Si at the scales \i = /xo and /x = We take 

= 4.9 GeV, m c = 1.5 GeV, and the light u or d constituent quark masses to be 0.3 GeV. Here 
L stands for longitudinal and T for transverse. 



Meson 


r 




\ z )n=Ho 




^z/2 








L 


T 


L 


T 




0.50 


4.5 


0.61 


0.62 


0.48 


0.48 


B* 


0.23 


7.9 


0.73 


0.73 


0.61 


0.61 


D* 


0.17 


2.1 


0.77 


0.77 


0.50 


0.50 


B* 


0.058 


5.5 


0.87 


0.87 


0.70 


0.70 



Table 2: A table showing where the curves D L {t) and D T (t) cross-over, the average (t) = (p±)/M, 
and the average (p±) for the b quark fragmenting into B* and B*, and for the c quark fragmenting 
into D* and J/ip. The mass of each meson is taken to be the sum of the constituent quark masses 
(m b = 4.9, m c = 1.5, m u>d = 0.3 GeV). 





r 


t (cross-over) 




(*> 


(P±> (GeV) 








L 


T 


L T 




0.5 


1.6 


1.8 


0.86 


5.5 2.6 


K 


0.23 


1.5 


1.1 


0.61 


7.0 3.9 


D* 


0.17 


1.6 


0.86 


0.52 


1.5 0.94 


B* 


0.058 


2.0 


0.41 


0.30 


2.1 1.5 
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